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Abstract. The Hermite functions are an orthonormalbasis of the space of 
square integrable functions with favourable approximation properties. Al- 
lowing for a flexible localization in position and momentum, the Hagedorn 
wavepackets generalize the Hermite functions also to several dimensions. Us- 
ing Hagedorn's raising and lowering operators, we derive explicit formulas and 
recurrence relations for the Wigner and FBI transform of the wavepackets and 
show their relation to the Laguerre polyomials. 



1. Introduction 

The Hermite functions are an important member of the family of special func- 
tions. They form an orthonormal basis of the space of square integrable functions 
on the real line. The Hermite functions are eigenfunctions of the Fourier transform 
and of the harmonic oscillator. They can be generated either by raising and lower- 
ing operators or by a three-term recurrence relation. And one could name more of 
their distinguished properties. 

The Hagedorn wavepackets |H81[ |H85 ( H98J generalize the Hermite functions 
to several space dimensions while adding more flexibiliy in terms of position and 
momentum localization. Moreover, they allow for a positive scale parameter e > 
setting the wavepackets' width at the order of y/e and their wavelength at the order 
of s. In particular, they are constructed from the complex Gaussian wavepacket 

<f%(x) = (ne)-^ det(g)- 1 / 2 exp(£(x - q) T PQ-\x - q) + |p T (x - q)) 

centered in position q G M. d and momentum p 6 K d , where the real and imaginary 
parts of the complex matrices Q,P E C dxd satisfy a symplecticity condition. In 
this way, the Gaussian has a position density |¥>o( x )l 2 > which is proportional to the 
multivariate normal distribution with mean q and covariance matrix %QQ* . The 
corresponding momentum distribution |J re tp e (^)| 2 , defined via the e-scaled Fourier 
transform 

= (2ne)- d/2 [ <p(x)e- lxT t/ £ dx, £ € K d , 



is proportional to the normal distribution with mean p and covariance |PP*. The 
familiar raising operator of the Hermite functions is generalized to 

At = yfc ( p * (*-«)- Q* (-«V a - p)) , 

and we obtain the fcth Hagedorn wavepacket by the fcfold application of the 
raising operator to the initial Gaussian, that is, 
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In view of a construction process, which balances the position and the momen- 
tum parameters, we ask for transforms of the Hagedorn wavepackets, which treat 
position and momentum or time and frequency variables simultaneously. That is, 
we aim at the Wigner transform and the Fourier-Bros-Iagolnitzer (FBI) transform 
of the Hagedorn wavepackets, paying due attention to the ladder operator approach. 

The e-scaled Wigner transform of two Schwartz functions <p, ijj : R d — >• C is the 
Fourier transform of their correlation function, 

W £ {^ iP){x, = (2ire)- d [ ?p(x + \)^{x ~ D^^dy, (x, e M 2d , 

see for example |Fol Chapter 1.8] for a discussion of basic properties or |Br] for the 
interpretation as a musical score. The Wigner transform maps Schwartz functions 
to Schwartz functions on phase space, while respecting orthogonality in the sense 
that 

(1) (W e (^i,^i),W e (</J2,-02))L2(R^) = (^l,^2)£2 (R rf)(^l,V>2)l,2(K<i) 

for all Schwartz functions (pi,if2,ipi,ip2 ■ R d — > C. On the diagonal, the Wigner 
function W 6 (y) = W £ (<p, if) is real- valued with the position and momentum density 
as marginals, 

\p(x)\ 2 =[ \^m?=l W £ (<p)(x,Odx. 

However, W £ (ip) is not a probability density on phase space, since it might attain 
negative values. For example, odd functions ip satisfy W £ (tp)(0, 0) = — (27re) _d ||(y9|| 2 , 
and also the Wigner transforms of the Hagedorn wavepackets, except for the initial 
Gaussian ip^, attain negative values. This lack of positivity can be cured by the 
convolution with the properly e-scaled Gaussian phase space function G e (z) — 
(Tre)~ d exp(— -j|z| 2 ), z £ M. 2d . The resulting positive transform 

H e (<p) = G £ *W £ {p) 

is the so-called Husimi transform H £ ((p) ■ M. 2d — > [0, oof. The Husimi transform 
can also be deduced from the Fourier-Bros-Iagolnitzer (FBI) transform, which is 
defined as the inner product with the Gaussian wavepacket 

9l ( (y) = (*e)- d/i exp(-^|y - x\ 2 + §£ T (y - x)) , ye R d , 

centered in the phase space point (x,£) £ R 2d . That is, 

T e (<p)(x,0 = {2-Ke)- d ' 2 ^,9l^), M e K M , 

see for example |Fol Chapter 3.3]. Then, the Husimi transform appears as the 
modulus squared of the FBI transform, that is, H £ (<p)(x,(;) = \T £ (ip)(x, £)| 2 for all 
(x,£) € M. 2d , which immediately reveals positivity. 

Our study of the Hagedorn wavepackets from the phase space point of view pro- 
ceeds as follows. In Section ^\ we give a brief account on Hermite functions and 
recall a polynomial sum rule for computing the Wigner and FBI transform of the 
Hermite functions. In Section we review the construction process of the Hage- 
dorn wavepackets, discuss their relation to the Hermite and Laguerre polyomials, 
and provide a Rodriguez type formula. Section SQ] then derives explicit formulas 
for the Wigner and FBI transform of the Hagedorn wavepackets and recasts the 
three-term recurrence relation of the Hagedorn wavepackets on the Wigner level. 
In the appendix, Section SjX] presents another characterization of the Hagedorn 
wavepackets based on the polar decomposition of one of the width matrices, while 
SjB] reformulates the ladder operators in Weyl quantization. 
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2. Hermite functions 

Hcrmitc functions are Hermite polynomials times a Gaussian. They can be 
generated from the Gaussian 

(2) tpo(x) = n- 1/4 cxp(-±x 2 ), i£l, 
using the ladder operators = -^{x — V^), 

^ fe +! = 7m Af(pk > k e N - 

The formal adjoint A = -^{x + \7 x ) of the ladder operator A^ allows to descend 
within the Hermite functions, 

(p k = -jj=fA(p k+1 , fceN. 

Remark 1. The Hermite functions are eigenf unctions of the Fourier transform 
and the harmonic oscillator ^(AA^ + A^A) — |(— A x + x 2 ), 

T\ k = (-*)**>*> h( AAl + ^A)tp k = (k+ \)<p k 

for all k £ N. Moreover, the family {ip k \ k € N} forms an orthonormal basis 
of the Hilbert space L 2 (JBL), see for example |T1[ §7.8]. They enjoy the following 
approximation property: Let K £ N, s < K and f : M. d — > C a Schwartz function. 
Then, 

||/- X>*,/>v*L 9(R) < (x(^-i)---(^- s + i))- 1/2 ||^7ll L2(R) , 

k<K 

see |Lul Theorem 1.2]. 

2.1. Hermite polynomials. Alternatively, the Hermite functions can be written 
as 

ip k (x) = —^=h k (x)ipo(x), x E K, 
V2 k k\ 

with 

(3) h k (x) = exp(x 2 )(~£) k exp(-x 2 ) 

the fcth Hermite polynomial. Starting from ho = 1, the Hermite polynomials can 
be generated by repeated application of the ladder operator = 2x — V x , 

h k+1 = B^h kl k e N, 

or from the three-term recurrence relation 

h k+ i{x) — 2xh k (x) - 2khk-i(x), k > 1. 
The orthonormality of the Hermite functions implies for the Hermite polynomials 

(5) f h k (x)h l (x)e~ x2 dx = Vn2 k k\5 u , fc,l € N. 

Jr 
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2.2. Integral formulas. The Hermite polynomials satisfy several beautiful inte- 
gral formulas. Those, which we employ for the phase space transformation of the 
Hermite functions, can be deduced from the following sum rule 

(6) h k (x + z)=Y / ( k )(^) k - j h j (x), x,ze£, 

j=o W 

which is due to [E] ■ We refer to Proposition [3] later on for a proof. 

Proposition 1. Let k < I and h k and hi be the kth and Ith Hermite polynomial. 
Then, for Z\,Zi G C, 

(7) f h k (x + z 1 )h l (x + z 2 )e- x2 dx = ^2 l k\z l 2 - k L<i- k \-2z 1 z 2 ), 

JR 

where 

j=0 V ' 



are the Laguerre polynomials associated with 7 G R. In particular, 
(9) / h t (x + z 2 )c~ x2 dx = y^2 l z 2 . 

JR 

Proof. We write 

hk(x + z\)hi{x + z 2 )e~ x dx 



j=0j'=0 

From the orthogonality condition Q we then deduce 

k 



J^h k (x + z 1 )h l (x + z 2 )e- x2 dx = Q ^(2zi) k - j (2z 2 ) l -^^j\ 



3=0 

k 



^^E ,,^,),^! 2 ^ = ^klz^Lt k \~2z lZ2 ). 



□ 



2.3. Phase space transforms. The Hermite-Laguerre connection of PropositionQ] 

translates to the Wigner transform of Hermite functions. For alternative proofs, 
see [12 Chapter 1.9] or [Tnl Chapter 1.3]. 



Corollary 1. If (p k ,<pi are the kth and the Ith Hermite function, then the Wigner 
function is 



(^y^y|#-fc e -M 2 4'- fc )( 2 | z | 2 ), k < I, 
t^^2^^ z k~i e -\ Z \ 2 L ( k - l ) {2 |z| 2 ), I < k, 
with z = x + i£ for x, £ G R. In particular, 

W^ k )(x,0= { -^e-^L{°\2\zf). 
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Proof. We compute 



where we have used hi(—x) = (—l) hi(x). Changing the variable as y/2 — i£ = 77, we 
perform a contour integration in the complex plane. Analyticity and exponential 
decay of the integrand then provide 

The integral formula (J7J| concludes the proof. □ 

The integral formula © of Proposition [1] provides the FBI and the Husimi 
transform of Hermite functions, see [PT] §2]. 

Corollary 2. Let ip^ be the kth Hermite function. Then, the FBI transform is 

V7r2 fc+1 fc! 

with z = x + i£. Consequently, the Husimi transform is 

Proof. We apply the integral formula |S| to obtain 

TH^ix^) = 1— f h k {y)e-fr e -hv-*)'e*<v-*)dy 
7rV2 fc+1 fc! Jm 

_i , 1 

e 2 x « e 4 



7rV2^+ T fc! 
_i , 1 
e z xi e 4 



*? a -z(* 2 +e 2 ) 



hk{y)e -(y-2^) dy 

(x + iO k - 



□ 



Remark 2. We no£e £/iat both the Wigner transform W 1 and the Husimi 
transform Tl l (ipk) of the kth Hermite function tpf. are radially symmetric, that is, 
they are functions of the energy variable \z\ 2 = x 2 + £ 2 , see also |HC[ §2.4]. 

3. HAGEDORN WAVEPACKETS 

In [H98J , George Hagedorn devised parametrized ladder operators which allow for 
a beautiful generalization of the Hermite functions also in several space dimensions: 
Let e > be a scale parameter, q,p g R d and C = C T G C dxd a complex symmetric 
matrix with positive definite imaginary part Im(C) > 0. 

Remark 3. Any complex symmetric matrix C € <C dxd with Im(C) > can be 
written as C — PQ^ 1 , where P,Q € C dxd are invertible matrices satisfying 

(10) Q T P-P T Q = 0, Q*P- P*Q = 2ild. 

This condition is equivalent to 
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Conversely, any pair of matrices P,Q G i£ dxd satisfying ilO\) is invertible and 
defines via C — PQ^ 1 a complex symmetric matrix with 

Im(C) = {QQ*)- 1 > 0, 

see |Lu[ §5, Lemma 1.1]. The applications of the Hagedorn wavepackets especially 
to quantum dynamics emphasize the importance of allowing both matrices Q and P 
to be complex. We note, that Hagedorn uses A and iB for Q and P in his work, 
while we adopt the more recent notation of |FGLl ILu] . 

Let Q,P G C dxd be matrices satisfying (]TU]) . Then, the Hermite ladder operators 
= -^{x — V a ) and A = -^{x + V x ) generalize as 

A%p,Q,P] = j=(P*(x-q)-Q*(-ieV x -p)), 

A[q,p,Q,P] = -^=(p T ( x - q )-Q T (-ieV x -p)). 

Before employing this ladder, we summarize some useful linear algebra. 

Lemma 1. Let e > 0, q,p G M rf , and Q,P G C rfxd (, e matrices satisfying U0\) . 
Then, QQ* and PP* are real symmetric matrices. Moreover, the components of 
A^[q,p,Q,P] = (Aj)j =1 commute, that is, 

4>4] = 44-44 = ° 

for j, k = l,...,d. 

Proof. Since C = PQ^ 1 is complex symmetric with Im(C) = (QQ*) -1 > 0, we 
have QQ* G R dxd and QQ* = (QQ*) T . Moreover, 

PP* = CQQ*C = Rc(C)Im(C)- 1 Re(C) + Im(C) G R dxd , PP* = (PP*) T . 



Finally, 



a a 

1=1 m=l 
d 

l,m=l 

-P*Q + Q*P)^ = 



2 V ■* 1 -e ~ Jjk 

due to the canonical commutator relation -^[xj , —ied Xk ] = Sjk and the matrix 
condition ([TU]) . 



□ 



The next step, is to generalize the zeroth order Hermite function ipo as the 
complex Gaussian wavepacket 

V>o( x ) = Vo[q>P>Qi p ]( x ) 

= (vre)-^ 4 det(Q) -1 / 2 exp(£(x - q) T PQ~\x - q) + \p T {x - q)) 

for x G R d . We note, that the vectors q,p G R d and the matrices QQ*,PP* £ R dxd 
provide the centers and the width of the corresponding position and momentum 
densities, 

\<pl(x)\ 2 = (^)^/ 2 det(Q)- 1 exp(-i(x-< Z ) T (QQ*) -1 ( 2 ;- g )), 

\F e <P e (0\ 2 = M- d/2 det(P)- 1 exp(-i(e-p) T (PP*)- 1 (e-rt) 

for x, £ G R d . With the ladder operator A^ — A^ [q, p,Q,P], the Hagedorn wavepack- 
ets tp e k — ipt [q,p, Q, P] are defined recursively via 



(11) 



Vk+e, 



=4^1, 



k G N d , 
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where ej £ N d denotes the jth unit vector for j = 1 , . . . , d. The formal adjoint of 
A$ allows to descend, 

The harmonic oscillator equation of the Hermite functions generalizes to 

d 

That is, every Hagedorn wavepacket ipf,, k £ N d , is an eigenf unction of an harmonic 
oscillator for the eigenvalue |fc| + \, see |H98[ Theorem 3.3]. 

Remark 4. The scaled Fourier transform of the kth Hagedorn wavepacket satisfies 

F s <p%[q,p,Q,P] = (-ipe- i * T <*/ e <p%\p,-q,P,-Q}, k £ N d , 

see |H98i §2] for an effortless ladder proof. 

The Hagedorn wave packets behave conveniently under modulations and trans- 
lations, which will be useful for computing their Wigner and FBI transform. 

Lemma 2. Let e > 0. Let ip,ip : M. d — >• C be Schwartz functions and 
for (x,£) £ R 2d . Then, 

w e &,ip)(x,o = w e (<p x ,t,ip x<s )(o,o), t s W(x,o = r e (^,c)(o,o). 

Moreover, for k £ N d , q,p £ R d , and Q,P £ C dxd satisfying \1(J\) we have 

(ip%[q,p,Q,P}) x ^ = e-^ T ^ip%[q-x,p-^Q,P]. 
Proof. We check 

W% x , 6 ^)(0,0) = (2ne)- d ( ^(|)^(-|)^ 

= (2ire)- d f Tp{x + %)xp(x - %)e iyT ^ e dy 

and 

T e (ip)(x,£) = {2Ks)- d,2 {T:e)- d/A f ^(y^-^-^+e^^-^dy 

JR d 

= (2ire)- d/2 (ire)- d/4 f ^{y + x)e-h^+^ Ty dy 

JM d 

= r £ (v^,c)(o,o). 

Moreover, (ifi £ [q,p, Q, P]) x>£ = e~^ T {q ~ x) ip £ [q -x,p-£,Q, P] and 

(A* [q,p, Q, P]tp) x£ = A^q~x,p-Z, Q, P]^ 
for all Schwartz functions cp. Therefore, for all k £ N d , 

{<p%[q,p,Q,P]) x ^ = -^{{A^[ q ,p,Q,P]) k ^[q,p,Q,P]) xi 

= e-^ T to-*)<fi[q-x,p-S,Q,P]. 

a 
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3.1. Hermite polynomials. The Hagedorn wavepackets are known to satisfy the 
three-term recurrence relation 

(12) (yj^+iip% +e .tx)) d ^ = 

'-'*-\x - qM(x) - Q-'Q (yk~w%- ej (x)) d j=i , 

see |Lul Chapter V.2]. Therefore, (p k = ip k [q,p, Q, P] is the product of a polynomial 
of degree |fc| with the complex Gaussian ip^. However, we can deduce some more 
information. 

Proposition 2. Let e > 0, q,p G M. d , and Q,P e C dxd satisfy fW^). The kth 
Hagedorn wavepacket cpf, = ip £ k [q,p,Q, P], k G N d , can be written as 

(13) <f%{x) = -^=p%{x)^l{xl x e R d , 

where p| is a multivariate polynomial of degree \ k\ generated by the recursion p^ — 1, 
V e k+e3 = B]p% for j = 1, . . . , d, with 

B t = -^Q-^ 0Pe (x - (?) - -tQ*(-ieV x ). 

The components of = (Bj) d =1 commute, that is, [Bj,Bt,] = for j,j' — 
1, . . . , d. Moreover, 

(14) Pl(-x) = {-l)Wfi(x + 2q), xeR d . 
IfQ = Q T e M. dxd , then 

d 

(15) Pl(x) = l[h kj (±(Q- 1 (x-q)) ] ), xeR d . 

3=1 

Proof. We denote A' = A'[q,p, Q, P]. We first consider the special case Q = Q T G 
M. dxd and compute 

d 

Q*(-ieV x -p) J| h ki {Vj) <p e ( x ) 

3=1 

= -i^!Q*Q- T : cp E (x) 

l^WIlj^ h kj {Vj)) 

d 

+Q*C(x-q)l[h kj (y j )tp s (x) 

3=1 

with C = PQ- 1 . Since P*Q - Q*CQ = -2ild, we obtain 

d ( (%i^i (yi) - h 'k t (vi)) Uj^i h k 3 (yj 



A^Y[h kj ( yj )^ (x) = 



3=1 



ipl{x) 



1 

v/2 



\{2y d h kd (y d ) - h' kd {y d )) Y\^ d h kj (y J ) / 

^itife/ilIlj/i^feT 

ipl(x). 



Assuming, that the claimed identity (|T5|) holds for k G N d , we derive 
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For general Q £ C dxd , we observe that 

Q*(-»W a -p)p%{x) V %{x) = {Q*(-ieV x )p%) (x)<p s (x) + Q*C(x - q)jf h (x)<pftx) 
and 

Aipl(x)tpl(x) = -^((P*-Q*C)(x-q)p%(x)-Q*(-ieV x )p%(x))<pl(x) 
= 72 {tS Q ~^ x - ~ j= e Q*(-ieV x )p%(x)) v g(x). 

Assuming that equation ([T31 holds for k £ N d , we conclude 

Moreover, 

d 

/ , rn — 1 

= 2e{Q- 1 Q-Q*Q- T ) rr = 0, 

where we have used that [xi, —ied Xm ] = ieSi m and QQ* — (QQ*) T — QQ T ■ For 
proving the symmetry relation, we argue once more inductively. We have 

Pf+e 3 .(-*) = ^{Q-\-x-q)),pl{-x)-^e{{Q*S/ x ) j p%){-x) 

= ^(Q-\- X - q )),(-l)\%l( x + 2q) 

+Vi(Q*V x )j(-l)^p%(x + 2q) 

= (-l)W +1 pl +ej (x + 2 q ). 

□ 

Remark 5. In the univariate case, Proposition^ and[7J provide the orginal defini- 
tion of tp £ k [q,p,Q,P] given in |H81l §1], 

V%[q,P,Q,P](x)= ( ^| ^R. 

for Q = Q T G R. dxrf , the description of Hagedorn wavepackets as rotated, scaled 
versions of harmonic oscillator eigenfunctions is mentioned in |H851 Remark 4 of 
§1]. For general Q £ C dxd , there is a more complicated relation, see Proposition^ 

3.2. Laguerre polynomials. The polynomial ladder of Proposition [2] allows to 
generalize the Hermite polynomial's sum rule ^ and the integral connection to the 
Laguerre polynomials. A similar argumentation for generalized univariate Hermite 
polynomials can be found in |DMI Proposition 2]. 

Proposition 3. Let e > 0, q,p £ R d , and Q,P £ C dxd satisfy {77$. Let k, I £ N d , 

k < I, and pt, p\ be the kth and lth polynomials defined in Proposition^ Then, 
for x,y,z £ C d , 

(i6) p%(x+z) = J2 (*) (^<r^) fe ~V» 

u<k ^ ' 

and 

~ d 

/ ¥ k {x + y)pt{x + z)\ v l{x)\ 2 dx = J] C kjh ( A. (Q=*v) j: MQ^z),) 
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with 

r ( , / 2™m!C n ~ m ^ m) (-2r/C), m < n, 

Proof. For z £ C d we use the translation {r z f)(x) = f(x + z) and observe 

r z oBt = (B^ + ^ s Q- 1 z)o Tg . 

Inductively, we obtain 



Tz o(B^ k =(B^ + ^Q^z) or, 



for all k € N d , which applied to P q = 1 yields the claimed sum rule. Therefore, 

,j<-i.,j'<-i \ / \ / 



min(fei,/i) min(fc <J ,/ d ) 



E - E (t -^'- y)M ' M (jb^r(jfeo-^ 

1/1=0 y d =0 v y v y 

where we have used the orthogonality relation 

If kj < Z j , then 

= fc,!2^ (^(Q-i^)^" fc3 '£g-^)(-2 ( ^Ty )i( Q-i^ 

by the monomial representation of the Laguerre polynomials ([5]). The analogous 
argument for kj > lj concludes the proof. □ 

3.3. Rodriguez formula. The Hagedorn wavepackets' three-term recurrence for- 
mula (|12p rewrites on the polynomial level as 

(17) (p% +e . (x))^ = ^Q~\x - q)pt(x) - 2Q- l Q (k jP e k _ e . (x)Y j=i . 

Together with the following Rodriguez-type formula we obtain another integration 
formula for the polynomials. 

Proposition 4. Let e > 0, q,p £ M. d , and Q,P G C dxd satisfy [TQ). Let k e N d , 

and p| be the kth polynomial defined in Proposition^ Then, 

p%(x) = \yl{x)\-\-V~sQ*V x ) k \vl{x)\ 2 , x e M d . 

Let M — \{iC + Id). Then, Id + Q T A1Q is invertible, and we have for all z G C d 

[ p%{x + z)e-^ T ^ +M ^dx = V ( k ) U0^z) k ^ cu 
Jm« „< fe \vj ^ ' 

with 

cu = 0, (c v+ej )j=i = -2(Id + Q T MQ)- l Q T MQ{v c v ^) d =1 
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for v £ N rf with \v\ odd. Moreover, M = is equivalent to C = ild. In this case, 
[ p%(x + z)e-^ x ~ q)Tlm(CKx - q) dx = {nef/ 2 det(Q) (^Q^zY . 

Proof. We set q%(x) = \(Pq(x)\~ 2 (— *fsQ*V x ) k \tpl(x)\ 2 and verify, that q e k satisfies 
the three-term recurrence (|f 7| . By the Leibniz rule, 



\vl(x)r 2 {-V-eQ*V x ) k ^Q-\x qW {x)\ 2 



koWl" 2 E ( ^ {-V~eQ*V x y^ E Q-\x q)(-V~sQ*S7 x ) k -"\^ ]2 



W Q {x)\- 2 ^Q-\x - q)(-V~eQ*V x ) k \^ (x)\ 2 
d 

2Q- 1 \^(x)\- 2 Y J k ] {{Q*^ x ) ] x) (-^Q*V x ) fc -^|^( 
i=i 

^ e Q~\x - q)ql(x) - 2Q- 1 Q (k jq %_ e . (z))^ . 



x)\ 2 



Moreover, 

(Id + Q T MQf = fid + \Q*CQ = fid + \Q*P = Id + §Re(Q*P), 

where we have used the matrix condition (fTu| . Hence, Id + Q T MQ is invertible. 
For proving the claimed integral formula, we use the sum rule (|16p and obtain 

I W k (x + z)e-^ x ^ T ^ +M ^dx = £ ( fc ) (^Q^z) fe ^C 

with 

<V = / W v {xKl (x ~ q)T(lm(C)+M){x - q) dx. 



By the symmetry relation Q3]), c„ = for odd. Moreover, by the Rodriguez 
formula and the three-term recurrence, 

(( _^gT Vxr+ e i) ^ i |^ (x)| 2^ 

-^Q T M(x - q)e^ {x ' q)TM{x - q \-^eQ T V x r\^{x)\ 2 dx 
-Q T MQ f -^Q^(x -qM(x)e-^ x - q)T ^ c)+M){x - q) dx 



-Q T MQ(cl +e] )* =1 - 2Q T MQ(^cU e ^ =1 . 



□ 



Remark 6. In the univariate case, the previous integral formula allows for simpli- 
fications. We set i = Im(C) + M and a 2 = Q 2 M/{1 + \Q\ 2 M) to obtain 



Jmx +z) e-^ x - q Ux= e (*) 



with 

Cj = -2a 2 (2i-l)c J _ 1 =2^-a 2 ) i (2i-l)-(2i-3)---l-co 



V^T(-a 2 ) J (2j)!/j. 
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The monomial representation of the Hermite polynomials then implies 

p%(x + z)e~~^ (x ~ q) dx = y/neaia^hk {-^Q^z^j . 
For the Hermite polynomials, this formula is due to |Fe| . 

4. Hagedorn wavepackets in phase space 

Our studies so far have provided two integral formulas for the polynomial part 
of the Hagedorn wavepackets. Now we apply them for computing the Wigner and 
FBI transform. 

4.1. Wigner transform. Proposition's[3]integral connection to the Laguerre poly- 
nomials allows us to write the Wigner function of the Hagedorn wavepackets in 
terms of Gaussians and Laguerre polynomials depending on the complex vector 

z(x, = -i {P T (x -q)- Q T (i -p)), x^e R d . 

Theorem 1. Let s > 0, q,p E R d and Q,P £ C dxd satisfy 1777)) . Then the scaled 
Wigner function of the kth and the Ith Hagedorn wavepacket ip^. = ip e k [q,p, Q, P] 
and (pf — (pf[q,p,Q, P], k,l £ N d , satisfies 

(is) = ^r^^mm. n*w**> 

with z=-i (P T (x -q)- Q T (£ - p)) for (x, £) e R 2d and 

f 2"m! C n - m i^" m) (2|C| 2 ) , m < n, 

mMU \ 2 m nl (-C) m ~"Ll m -" ) (2|C| 2 ), n<m. 



In particular, 



j'=i 

Proof. We first study the Wigner transform in the origin. By Proposition [21 

w £ (^f)(o,o) 

V2l' c l+l i lfc!Z! ii-i 
_. (-i)IH /■ _ _ 

= ( 7rE ) /f ,|,M|,|, ,„ / pUv - q)Pi(v + q)<Po(v - q)<Po(-v + v)dv- 

V2l fc l+I'lfc!i! JR<i 

The idea is to simplify the integrand such that the Laguerre connection of Propo- 
sition [3] can be applied. We compute 

fo(v~ q)<pI (-v + q) 

= (ne)- d ^ 2 det(Q)- 1 e-i( I '- 2 «) T ^'"- 2 ^+i' )Tc "-f f T (^«) 

= ( 7r£ )-d/2 det (Q)-l e -7') T Ini(C) I ?+fr ; T Im(C) ? -fg T Im(C)9+f ( V -q) T (Be(C)q-p) _ 

The complex vector 

z=-i (pT(-q)-QT { - p) ) 

satisfies z — iQ T (Cq — p), since C is complex symmetric. Due to Im(C) = {QQ*)~ l 
and QQ* = (QQ*) T = QQ T , we have Qz = -q + i\Q\ 2 (Re(C)q - p). Moreover, 

\z\ 2 = z T z = q T lm(C)q + (Re(C)q - pf\Q\ 2 (Re{C)q - p). 
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We set w = —ilm(Qz) = — i\Q 2 \(Rc(C)q — p) and compute 

{ip^ + w)\ 2 = (7re)- d / 2 det(Q)- 1 e-^" +t "-« )Tlm(c ' )( " + " , - <?) 

= f n£ \-d/2 det( -Q)-l e -ig T Im(C)g+i(Rc(C)g-p) r |Q| 2 (Rc(C)g-j 9 ) 



Therefore, 



e -i t? r Im(C)r,+ f J ? r Im(C)g+f (, 1 ~q) T (Rc(C)q-p) _ 



<p e (v - q)<pU-v + q) =e ei 2| Vo(^ + HI 2 



and 

W £ (^,^)(0,0) 

= ^ ; n iii x i,i,,,r / P%(y-'w-q)p e l (y-w + q)\(p e (y)\ dy 

due to analyticity and exponential decay of the integrand. By Proposition [51 

W%1,^)(0,0) 

= (7T £ )- d ( ~ 1) '" e = IT ( l Zj ,-lzj) , 

where we have used that —q — w = Qz and </ — u> = — Qz. For arbitrary (x, £) G R 2d , 
we use Lemma [5] and obtain 

W%|, <pf = [<z -x,p-£,Q, P],tf{q -x,p-t,Q, P])(0, 0). 

□ 

To generalize the observation, that the Hermite function's Wigner function only 
depends on the energy variable \z\ 2 , see Remark^ one has to combine the Hagedorn 
wavepackets for the |fc|th eigenspace. We set $f fe , = $f fc , [q,p, Q, P], 

*f*l (*) = K (*), ■ • • . vf»)/ l+ -+/ d =i*l e x e K d , 

and move to the associated matrix- valued Wigner transform 

Corollary 3. Let e > 0, q,p £ R d and Q,P E C dxd satisfy {37$. Then, the 
Hagedorn wavepackets' vector for the \k\th eigenspace, $^ fc | = 3>f fe i [q,p, Q, P], k G 
N d , satisfies 

with z=-i {P T {x -q)- Q T (£, - p)) for (x, f ) G R 2d . 
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Proof. 

tr(w £ ($f fe| ,$f fc| )(x,£) 



(_l)l*l 



|z|2 E 4 0) (lM 2 )--^L 0) (iM 2 ) 

h+—+ld=\k\ 



(-1)1*1 l lz{ 2 (d-l), 2| , 2) 

by the Laguerre polynomials' summation theorem |Ra[ §119]. □ 



4.2. Three-term recurrence relation. The Hermite functions' ladder operators 
can be translated to the Wigner function level, see |Tnl Theorem 1.3.3]. The same is 
possible for the Hagedorn wavepackets and provides a useful three-term recurrence 
relation in phase space. 

Theorem 2. Let e > 0, q,p € R d , and Q,P E C dxd satisfy (W$). Then, the 
Wigner transform WL — VV e (ip|,yf) of the kth and Ith Hagedorn wavepackets 
tp% = V>\ fe' P> Q> P\ an d <Pi — <Pi P>Q, P], k,l e N d , satisfies 



for j = 1, . . . , d with 

K ^ = ^(^ T (2(x- g ) + H £ v c ))-o T (2(e-p)-H £ v,))), 

K = --^_(P*(2(x- ? ) + (-feV c ))-Q*(2(C-p)-(-feV B ))) 

and 

Lt = ^(P*(2(x-q)-(-ieV i ))-Q*(2^-p) + (-ieV x ))), 
L = -^(P T Mx-<l)-(-^ ( ))-Q T (2(t-p) + (-zeV x ))). 
Proof. We compute 



1 tT. 



(-ieV x )W s kl (x,0 = -(27rs)- d (-ieV x -p)<f%(x+$)<fi(x-*)e^ *dy 



+ (2ne)- d I ?i{x + |)(-*eV x - p)<pf(x - fte^dy 



and 



(-zeV ? )>V fc e z (x,0 = (2ire)- d I (x + f - q)^ k {x + § M(x - § )e^^dy 



Since 

At = ^(^-ffJ-O'C-feV.-p)), 
A = --j=(p T ( x -q)-Q T {-isV x -p)), 

we obtain 

(p*c-fev c ) + q*{-iev x )) m&o = 

(2 TC )- d / (Aipl{x+ZM{x D + ^Kx+I^Vf^ - I)) eS^dj/ 

./Bti V / 



•'2? 
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and 

i 

V2e 



(P T H £ V ? ) + Q T (- !£ V X )) M&foO = 

./lttd V / 



Using A]ip% = y/kj + l<p% +e . and A 3 ip% = y/k~ip%_ e . , we arrive at 



Moreover, 

2(x-q)W s kl (x,Z) = (2ne)- d [ (x + § - q)^ k (x + *)<p?( x - ^e^dy 



Since £ee 4 y — — ieVj,ee 4 v , an integration by parts yields 

2(t-p)WZt(x,t) = {2*Ey d ( (-ieV x -p)tp%(x + %)tf(x-%) 



(2ire)- d / <pl(x+%)(-ieV x -p)<p!(x-%)est y dy, 



Consequently, 

^ {2P* (x-q)- 2Q* (£ - p)) mi = 

(2ne)- d f (-A<p%(x+%)tf(x k)+^+Y) A ^ l (x - f )) eU^dy 

.lad V / 



and 



jL (2P T (x - g) - 2Q T (£ - p)) Vyg, (x, 



(2^)- d / (-it^ + fM( a; -!) + ^(x+ f))et« T ^y 
This implies 



/2s 



(2P* (x- q )- 2Q* (e - P )) 3 m = - v^m- ej , t + Vh+^m^ , 
-^( 2pT ( x 'i)~ 2 Q T (^~p)),mi = -Vk^m+e^i+Vhm,^- 

Adding and subtracting the above raising and lowering identities for VVjy gives the 
operators K', K. U , and L. □ 

The phase space ladder allows to reformulate the Hagedorn wavepackets' three- 
term recurrence relation (|12p for the Wigner transform. 

Corollary 4. Let e > 0, q,p € M. d , and Q,P e i^dxd ^ e ma i r { ces satisfying kl(J\) . 
Then, the Wigner transform W|; = W s ((p kl (pf) of the kth and the Ith Hagedorn 
wavepackets tp k — (p k [q,p,Q, P] and cpf — (p £ [q,p,Q,P], k,l G N d , satisfies 

(Vk7+im + e 3 Ax,0) d j = 1 = + (v^WM-e>,0)" =1 

d i— / s d 



( v / ^+T>vi,i + e,(x,e)) i=l = sjizmiM + (V^m-ejA^o) , =1 

tuitfi z = -i (P T (a; - g) - Q T (£ - p)) /or (x, f ) € M 2rf . 
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Proof. We observe that 

(K]-L 3 )mi = ^{2P T {x-q)-2Q T {^p)) j Wt 



(L\-K 3 )W S U = ^L{2P*{x-q)-2Q*(£-p)) j W s H 

= 2^+lW e ktl+e . - 2yfcm-e j ,l- 



□ 



The numerical computation of the Wigner transform of a Schwartz function 
if; : R d — > C for a set of phase space points Cii • • • i Cw G is notoriously difficult, 
since a direct approach poses the numerical quadrature of N Fourier integrals in 
possibly high dimensions, 

W e (V>, 1>)(x, = {2^ey d [ + § )%l>{x - %)e iyT t/ e dy 



for (x,£) € {Cii ■ • ■ i Cjv}- The three-term recurrence of Corollary 2] might provide 
an efficient alternative, especially when a large number N 3> 1 of evaluations is 
required: Choosing suitable q,p € R d , Q,P e C dxd satisfying (fID|). ifeff, and the 
hyperbolic multi-index set 

d 

/C = {fceN d | n( 1 + fc j) ^ R }- 

3=1 

Then, one approximates 

^ ~ c fc <^|, c fc = / <p%(y)ip(y)dy 
fceJC ^ 

with ip € k — (p%[q,p, Q, P] for k € JC. With literally the same proof as for the Hermitc 
functions |Lu[ Theorem 1.5], the approximation error can be estimated as 

11^ - ^2 c kfl\\ L 2 {Rd) < C M X~ s/2 max || A[q,p, Q, P] a il>\\ L ^{^) 

for fixed s € N and every Schwartz function i\) : Mr — > C. From the Wigner 
transform's bilinearity and the orthogonality relation fT]) we deduce 

||W £ (^)-W £ (^ 2 )|| i2(R2t!) 

< ||W £ ((pi,(^l - <y9 2 )|| L 2 (R 2 t i ) + |1 W £ (^l ~ ¥?2, ^2)||i 2 (R 2d ) 

= \\fl - ^2|U 2 (R d ) (||< ( 5i|| L 2( Rc i) + \\ip 2 \\L 2 (R d )) 

for all Schwartz functions ipi,(f2 ■ R d — > C. Therefore, evaluating the Wigner 
function of ip via 



kJEK. 



together with the three-term recurrence relation of Corollary [?J we inherit the 
approximation accuracy of 0(K~ S / 2 ). 



4.3. FBI transform. Our general findings for the FBI transform of the Hagedorn 
wavepackets are less beautiful than those for the Wigner function. 
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Proposition 5. Let e > 0, q,p € R d and Q,P e C dxd satisfy (ffl\). Let tp e k = 

tp%[q,p,Q, P] be the kth Hagedorn wavepacket. Then the scaled FBI transform is 



(ire) d V2\ k \+ d kl 



l/<fc 



twi/i w = (iC + Id) 1 ((x — q) + i(£ — p)) for (x. £) g ]R 2d . In particular, if C = ild, 
then 



2A- 



(7re) d /V2l fc l+ d A;! 

and 

= (7r£)d2 | fc | +dfc , e-^ N2 

with z=(x-q)+i(£- p) for (x, £) € R 2d . 
Proof. We start for = (0,0) and obtain 

(7re) d V2l fc l+ d fc! •/jt" 

We compute 

-Uv ~ v) T c(y - i) - yty - 1) - Uv\ 2 

= -Teiv- q) T (iC + Id)(y -q)- ±(ip + q) T (y - q) - ±\q\ 2 
= -Uv-1- w) T {iC + ld){y -q-w) + ±w T (iC + ld)w - i|<?| 2 
with w = —{iC + Id) _1 (ip + q). Therefore, 

T%!)(0,0) = det( Q)- 1 / 2 ^ }q f + i. w T {{ c +u)w 



(7re) rf V2l fe l+ d fc! 

p|(y + w)e-^ {y - q)T{iU+w)(y - q) dy. 



Since ^(iC + Id) = Im(C) + M, we have by Proposition 0] 

r%f)(o,o) = 

det(Q)-V 2 e _^ + ^ wT{iU+u)wy (k\ / 
(Tre)V2l fc l+ d fc! ^ k ^^ 1 

For arbitrary (a;,£) G M 2d , we use Lemma [5] and obtain 

T%|)(x, = el^^T^f [g -x,p-£,Q, P])(0, 0). 



In the special case C = «Id, we have QQ* = Id, P = 0, id = \{{x — q) + i(£ — p)) 
and c v — for ^ 0. Consequently, 



with z = (x — 9) + z(£ — p). □ 
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Appendix A. Polar decomposition 



The following observations are not needed for computing the phase space trans- 
forms of the Hagedorn wavepackets. However, they shed further light on the relation 
to the Hermite polynomials in the general case of invertible Q G C dxd . 

Lemma 3. Let e > 0, q,p G K^, and Q,P G <C dxd be matrices satisfying U0\) . If 
we decompose Q = \Q\U* with \Q\ = (QQ*) 1 ^ 2 and U G C dxd a unitary matrix, 
then PQ~ l = {PU)\Q\~ l , and the pair \Q\, PU satisfies condition \1(A) . Moreover, 



(19) 

Proof. We observe 



A^[q,p,Q,P} = UA^[q,p,\Q\,PU}. 



= Q T P-P T Q = (\Q\U*) T P - P T \Q\U* 

= U- T (\Q\ T (PU)-(PUf\Q\)U*, 

2ild = Q*P- P*Q = (\Q\U*)*P- P*\Q\U* 

= U{\Q\*{PU)-{PU)*\Q\)U*. 

Therefore, \Q\ T (PU) - {PU) T \Q\ = and \Q\*(PU) - (PU)*\Q\ = 2ild. For the 
ladder operators we compute 

A%p, \Q\,PU] = -j=((PU)*o Ve (x - q) - |Q|*op E (£ - p)) 

= 7Te U ^ P * ^ X - «) - Q*°Pe(Z ~P)) = U*A^[q,p, Q, P]. 

□ 

The unitary relation (flU)) of the ladder operators provides information on the 
Hagedorn wavepackets associated with indizes of the same modulus: We first 
enumerate the multi-indizes k G N d of equal modulus redundantly by setting 

/>o = (0, . . . , 0) e N d , 

/ ^|fc|,i 



"\k\+i 



vec 



C'l 



v\k\,i + ed 



\k\ G N, 



V 



\ v \k\,d\>>\ + e i ' ' ' v \k\,dW + e dj 

such that t'lfci+i is a vector of length d' k ' +1 , whose entries are multi-indizes in N d . 
Then, we mark repeated occurrences of multi-inidizes by setting 

oo 3f < j : v\ k \ tj , = i>\ k \j, 



V \k\,3 



otherwise, 



for j = 1, . . . , This redundant book-keeping allows to reformulate the creation 
process (jilt on the level of the \k\i\\ eigenspace as 0q = (fg, 



V|*|+i 



vec 



/ 1 /it ,7tS 



We note that the normalization with 1 / \J{v\k\ j)i + 1 produces a zero whenever a 
multi-index is repeated in the enumeration. 

Proposition 6. Let e > 0, q, p G R d , and Q,P G C dxd satisfy We decompose 



= \Q\U* with unitary U G C 



dxd 



Then, 



\QV e 



T>dxd 



(20) 



<Pf k \ [q,P, Q, P](x) = U^^ k{ [q,p, \Q\,PU], 



where £/®l fe l = ?7® • • • ® Z7 G C d ' M xd ' M denotes the \k\-fold Kronecker product of U 
with itself. 



HAGEDORN WAVEPACKETS IN TIME-FREQUENCY AND PHASE SPACE 



19 



Proof. We set 

= A%p,Q,P], tf k{ =tf kl [q,p,Q,P], 
Dt = A%p,\Q\,PU], tf k \=tf k \[q,p,\Q\,PU]. 

By the relation 1(19)1. A] = (UD^)j = uj L> t with m,... ,u d £ C d the row vectors 
of 17. Moreover, for arbitrary W € c mxm , w € C m , and a e C d 

/ (u^a)Ww\ / (unai)V^u; + • • • + (ui d a d )Ww^ 



\(udiai)Ww H 1- (u dd a d )Ww 

fuixW ■■■ u ld W 



\(u T d a)Wwj 



\u dl W ■ ■ ■ u dd W j 
Assuming that the claimed identity ((20)) holds for |fc|, we therefore obtain 









f aiw\ 


1 




= U ®W 






\a d w) 




\a d w/ 



<p\k\+l 



( -7T= L TTT D i fc\ i \ 



^ (|fc|+1) V4, + i 



□ 



Appendix B. Weyl quantization 

The raising and lowering operators A^ and A as well as the generalized harmonic 
oscillator | (A T A^ + (A^) T A) of the Hagedorn wavepackets can be viewed as Weyl 
quantized operators obtained from smooth phase space functions of subquadratic 
growth a : R 2d ->• C d , 

(op £ (aM(x) = (27Te)- d f a(\(x + y), ^e^"^^)^ 



for Schwartz functions <p : M. d — > C. The emerging phase space function 

z(x,$ = -i (P T {x -q)- Q T (£,-p)) 

generalizes the complex number z(x, £) = x + i£, which characterizes the Wigner 
and FBI transform of the Hermite functions. It also appears in the Wigner function 
of the Hagedorn wavepackets. 

Lemma 4. Let e > 0, q,p £ M. d , and Q,P £ C dxd be matrices satisfying U0\) . Let 
A* = A%p,Q,P] and A = A[q,p,Q,P). Then, 

d 

^ = 7Te°^> A = 7%°^' 5 DM + 4 A i) = Te°Pe(\z\ 2 ) 

3=1 

with z(x,£) = -i (P T (x - (?) - Q T {£-p)) for x,£ £ R d and \z\ 2 = z T z. 
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Proof. We observe 

\z\ 2 = (x- qf\P\ 2 (x -q) + {£- p) T \Q\ 2 (t -p)-(x- q) T M{i - p) 
with M = PQ* +PQ T e R dxd . Moreover, 
\ (A T A^ + (A^ fA) = 

± (o Pe (x - q) T \P\ 2 op £ (x -q) + op £ (£ - pf\Q\ 2 op E (t - p) 
-o V£ (x-q) T Mo V£ (t-p) - {Mov £ (£-p)) T op £ (x-q)) 

and 

op £ (x - q) T Mop £ (£ - p) = op £ ((x-q) T M^-p))-e^tTM 7 
(Mop £ {Z-p)) T op £ {x-q) = op £ ((x^q) T M^-p))+e^tTM. 

□ 
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